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Abstract 

We generalize the Cottingham formula at finite (T 7^ 0) temperature by using 
the imaginary time formalism. The Cottingham formula gives the theoretical 
framework to compute the electromagnetic mass differences of the hadrons 
using a dispersion relation approach. It can be also used in other contexts, such 
as non leptonic weak decays, and its generalization to finite temperature might 
be useful in evaluating thermal effects in these processes. As an application 
we compute the vr"*" — ir^ mass difference at T 7^ 0; at small T we reproduce 
the behaviour found by other authors: 5m? {T) = 6m? {0) + 0{aT'^), while for 
moderate T, near the deconfinement temperature, we observe deviations from 
this behaviour. 
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1 Introduction 



The aim of this paper is to present the generahzation of the Cottingham 
formula to finite temperature (for a review of field theory at finite tem- 
perature and density see 0). The Cottingham formula allows the evaluation 
of the time-ordered product of two hadronic currents between hadronic states 
by a Dispersion Relation (DR) and its main application is in the evaluation 
of electromagnetic mass differences of hadrons. The use of the Cottingham 
formula in this context predates Quantum Chromo Dynamics (QCD) 
by the advent of QCD some modifications of the original formalism were 
adopted: for example, the problem of the convergence was settled by the 
use of an ultraviolet cut-off fi^ related to the renormalization procedure [Q. 
Furthermore, the use of effective field theories: chiral perturbation theory |^ 
and heavy quark effective theory [^], has introduced well defined theoretical 
schemes to evaluate the relevant Feynman diagrams. 

Besides the use in the context of the electromagnetic mass differences 
of hadrons, the Cottingham formula has ben also applied in other fields, 
mainly related to non leptonic weak processes: kaon decays [0, 0, — 
mixing p[, parity violations in nucleon interactions ||10[, -B-meson processes 
| |1 1|| ; therefore its extension to finite temperature may be of some interest to 
analyze the role of the thermal effects in these physical situations. Among 
the possible applications we shall pick up in this paper the 7r+ — 7r° mass 
difference at finite temperature, a subject that has received a continuous 



attention both in the past and in more recent times 13 



2 Cottingham formula at finite temperature 

As already stressed in the introduction, the main application of the Cotting- 
ham formula is for the calculation of the electromagnetic mass splitting of 
mesons. The mass shift of the meson M due to the electromagnetic interac- 
tion can be obtained by computing: 

2 J (27r)4g2 + ,e ^^'^^ ' ^ ' 

where the hadronic tensor 

T^"'{q,p) = 1 I (fx e-'''^M{p)\T{J>'{x)r{0)\M{p)) (2) 



describes the Compton scattering of a virtual photon of four momentum 
off the meson M of momentum p^ and J'^ is the electromagnetic current. The 



1 



Compton amplitude can be decomposed in terms of gauge invariant tensors 
as follows: 



where 



-9tiu + 




(3) 



(4) 



v = pq . 

The Lorentz invariant structure functions Ti(g^, v) and T2(g^, z/) depend, 

in the meson rest frame, by \q\ and q^ = — , where m is the meson mass. 

m 

For fixed the singularities in the complex plane are placed just below 
the positive real axis and just above the negative real axis; therefore one can 
perform a Wick rotation to the imaginary axis g° = ifc", without encounter- 
ing any singularity. After this transformation, the integration involves only 
space-like momenta for the photon: 



(5) 



After the change of variables (Bl) one obtains the Cottingham formula |]T| : 



5m 



IQn-^ Jo 



Q 



X 



-3Ti{-Q\imk^ 



'l'^]T2{-Q^^mk') 



(6) 



)2 

' max 



As discussed in and 0, the integration range is cut-off at Ql 
n"^. n is the Quantum Chromo Dynamics (QCD) renormalization mass scale, 
at which both the strong coupling constant and the quark masses nig have 
to be specified and roughly represents the onset of the scaling behaviour of 
QCD. As is well known the renormalization procedure introduces countert- 
erms which cancel the infinite contribution induced by virtual particles with 
momenta larger than /i: therefore its net effect is analogous to a cut-off of 
the integral at fi^. For approximate calculations there is a residue smooth 
dependence on fi, which in a complete calculation is exactly canceled by the 
/X— dependence of the renormalized quark masses and strong coupling con- 
stant. Typical values of fi are in the range of 1-3 GeV, corresponding to the 
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onset of the scaling behaviour of QCD and to a mass scale significantly larger 
than all the hadronic masses. 

We wish now to generalize this formula at finite temperature; we choose 
to work in the imaginary time formalism, which corresponds to substitute the 
integral over the energies with a discrete sum over the so called Matsubara 
frequencies Un = 2 7inT: 

This substitution corresponds to insert in (||) the factor 

+00 

27iT Y ^(^° - ^n) . (8) 

71= — 00 

We therefore obtain, from flBD and (E): 



~)2 _ ,,,2 



2 ■ 



where 

a = — 7^ 
2-nT 

[a] represents the maximum integer contained in the real number a. 

Eq.(|^) generalizes the Cottingham formula at T 7^ and is the starting 
point of the application to the vr"*" — 7r° mass difference at finite temperature 
to be discussed in the following section. 



(9) 



(10) 



3 The pion electromagnetic mass difference 
at finite temperature 

The invariant amplitudes Ti and T2 satisfy dispersion relations (DR) in the v 
variable. The DR for T2 is unsubtracted, while Ti requires one subtraction^. 



T,{q\u) = Ti g^O + — / -j^ ,2 ^ \ \ H 

IT Jo V v' ^ — 
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After the change of variables (|^) these equations become: 

^ du'^ Wii 



To 



Ti(-Q^O) 



m 



21.2 
fvn 



-Q 



2 U') 



,/2 



_|_ j^2 ^2 



+ 00 



/2 



i^2(-g 



2 zy') 



/2 



1,2). 



.(13) 
(14) 
(15) 



In order to evaluate the DR in the case of the pion electromagnetic mass 
difference we consider the contribution of the Born term (the vr meson itself) 
and the = 1~ resonances u. Other contributions might be in principle 
sizeable, however it is well known, since the work by Harari that the 
DR are sufficiently well convergent at T = and the first two polar terms 
in (|lT]) represent by far the largest contribution to the electromagnetic mass 
difference. This stems from the fact that the tt"*" — tt^ mass difference takes 
contribution from the Isospin= 2 mass term, and the time ordered prod- 
uct of the two electromagnetic currents in this case has no singularities for 
X — ^ (or, equivalently, for — s> oo), differently from other hadronic mass 
differences, such as the proton-neutron mass difference, that present a 1/x^ 
light-cone singularity. For the same reason the limit n ^ oo could also be 
taken in this case. 

To compute the different contributions to (|TT|) we consider the following 
matrix elements {q = p' — p): 



< n 



I 1^ 



TT 



~{p) > 



\JL 



F{q')ip + pr (16) 
7r\p)> = ih{q')e^'^'"^elq,p, (17) 

where e\ is the uj polarization vector and F, h are electromagnetic form 
factors. They can be written as follows: 

1 . . 



q'^/m 
h{0) 
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(19) 



1 — q^ /m\ 

my is a mass parameter that we identify with the p mass. Indeed both form 
factors can be obtained by assuming p dominance. This hypothesis allows to 
extrapolate the behaviour at finite temperature. Indeed at T 7^ 0, but small, 
one has 

,2 ■ 



F{q\T) 



2T2 
1 + — (?o 



X 



9pmT 



1 



q^ /rn? 




(20) 
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This expression takes into account the effective charge of a pion in thermal 
medium containing other pions at the equihbrium. Here = 93 MeV, go{x) 
is given by 



27r2 



dy 



y^x'^ + y'^ exp ( v^x^ + y'^) 



(21) 



while QpTTniT) and (7p^(T) are coupling constants for the pmr and p— photon 
vertex respectively, computed taking into account the effects of the thermal 
bath. They are obtained in the soft pion limit, which is sufficiently accurate 
for our purposes |]I4[. These effects (and the related /7r(T) dependence) 
have been computed by several authors and with different methods (see, for 
example []I1[,|TH] and references therein). The results are: 



(T) 



9pn 



9p-f 




The values of these constants at T = are related by the formula 

QpTTTrGp^t H 



(22) 
(23) 

(24) 



a relation that implements the idea of p dominance of the electromagnetic 
form factor. It is worth stressing that, independently of the thermal effects 
induced by electromagnetism, the pion mass gets a thermal self energy also 
by the strong interactions with the dilute pion gas. In general we shall not 
consider these effects since they are identical for vr"*" and 7r° and cancel in the 
difference. In one case the thermal self energy acts however as an infrared 
regulator (see below). We quote therefore this effect ||T3| : 



m^+,o(T) = 1 



(25) 



We also observe that, because of the interaction with the thermal bath, 
pions can couple directly to the photon without the intermediate virtual p 
state: this effect is at the origin of the last term in (^0]). This term, which 
only exists at T 7^ 0, gives a non-vanishing contribution to the form factor 
for — s> +00, which might be an artifact of the effective theory employed 
to compute it. However, as the integration is cut-off at p"^ and the results 
depends smoothly on p, there is no need to correct for this effect and we 
shall assume (120) as it stands. 
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Let us now turn to the form factor relative to the u intermediate state. 
We put h{0) = 2.6 GeV~"^, which can be derived by ^ 717 decay rate. As 
it will be clear in the sequel, a more precise determination of /i(g^) is not 
necessary because of the smaller role played by the u intermediate state in 
the calculation. For the same reason we omit to include thermal effects in 
the vertex wK'y that, to our knowledge, have not yet been computed. 

Using the matrix elements and the coupling constants already introduced, 
we can calculate the electromagnetic pion mass splitting. To start with, we 
consider the T = case. The contributions of the different terms to the DR 
follows; the subtraction term Ti(g^,0) is given by: 



T,{q\0) = -2F\q') + 



m?q^h?{q 



The two structure functions 1^1,2(9^5 ^) that appear in (|TB|) and 
j;iven by: 



W,{q',u) = +un{m'q'-ul)h\q')6{ 
W2{q\u) = 



2 2 



= -2m\^F\q^)5[v''-^ 
+ q'm'vnh\q')5(u'-ul] 



(26) 
5) are 

(27) 
(28) 



2 22 
m — m — q 

where vr = — . It is worth stressing that the Born term con- 
tributes only to W2{q'^, y), its contribution to Ti only being through Ti(g^, 0). 

These expressions can be used to compute the electromagnetic pion mass 
difference at T = 0. The results are reported in Table I, for /i = 2 and 
/i = 3 GeV and for fi = 00. It may be useful to stress that the numerical 
results are remarkably insensitive to variations of the cut-off. They show 
show that the pion contribution represents by far the dominant part of the 
mass difference. The small difference (around 10%) between the data and the 
theoretical result should be attributed to a few other poles in the dispersion 
relation, most notably the Oi resonance. 



Contribution 


/i= 2 GeV 


/i= 3 GeV 


fi = 00 


71 


3.56 


3.82 


4.05 




0.09 


0.11 


0.13 


Total 


3.65 


3.93 


4.18 



Table I. Contributions to 6m in MeV; dm^xp 



4.6 MeV. 
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In the chiral limit = we get, for /i = oo: 

= ^ , 29 

which gives, numerically, 6m = 711^^+ — m^o = 3.8 MeV, to be compared to 
the the experimental result Sme^p = 4.6 MeV. We can also compare it to the 



the result of Das et al. [0: 



^"^=^^lnh^ (30) 



that is obtained using the chiral limit and the Weinberg sum rules. The 
extra factor in (0) as compared to ( p9D is given by In £2 ) ^"^^ 

J TT J p J TV 

f 

numerically equal to 1.24 for — ^ = 1.66. 

Let us now consider the analogous calculation at T 7^ 0. We shall use in 



Ti{-Q',tmun) = -2F\-Q\T) 

_ ,„,„_«.)[^«^_^M±^] (31) 

T2{-Q ,imuJri) = Sm''F\-Q\T)- 



- m^h\-Q') , ^ , (32) 



where m?{T) is given in (^). By inserting these results in eq.(^ we can 
compute numerically the sums over the discrete energy. The numerical results 
will be discussed in the next Section. 



4 Electromagnetic mass difference in the chi- 
ral limit 

In the chiral limit the expression for the electromagnetic pion mass difference 
looks remarkably simple: 

60' /-m' dO'^ I 

^"^'(^) = S/o %F'^-Q'^T)2.T E v/Q^. (33) 

n = — [ml 



The behaviour of 5m? {T) at small T is as follows: 



(34) 



where drr? = 5rr?{T = 0) and A is a coefficient. 

The absence of the term linear in T can be proved explicitly by performing 
an asymptotic expansion in 1/T of 5m? {T). The term independent of T is 
given by 



limW(T) = ^/ dQ'F\-Q')\\ny- \/l - - 



47]-2 

6a 

3a /-A*' 



= £ I" dQ^F\-Q-) , (35) 

which coincides with 5m^, i.e. with the result obtained by putting directly 
T = from the very beginning. As to the coefficient of the term in T^, let 
us write it as follows: 

A = £1 + £2 , (36) 

as it arises from two different sources. The first term, £1, is obtained by 
expanding F'^{—Q'^,T) in T and taking into account that go{0) = 1/12: 



from this equation one gets the following contribution to 5m^{T) 

T2 



W(o) 



1 - 



3/i2 



m„ 



and, therefore. 



8/2 ^ ' V 3/i2 
The second term, i.e. C2, is obtained by putting 

F\-Q\T) = F\-Q^ 

and computing 



mt 



Co = lim 



6a r^'^ 
T-^o A^?f^ Jo 



rl "^^ ' r?2 ^ 



(37) 



(3^ 



(39) 



(40) 



(41) 



where, as before, a 



2ttT 



We have been unable to evaluate analytically 



this limit; however, when calculated numerically, for /i > 0.5 — 1 GeV, the 
limit is basically independent of ji and within a few percent is given by: 



£2 



ttq; 



Putting the two contributions together we get 



5m^{T) ^ 5m^{Q) 



6/2 V 4 



3 _^ + /i^ ^ ml + fi'^ 



(42) 



(43) 



holds 



We observe however that the behaviour of the last term in eq. ( 
only at small T (T < 100 MeV). For T > — , its T-dependence would be 
linear: 



3aT 



TT 



(44) 



but, given the range of the possible values of /i, such a behaviour cannot 
be reached, as the deconfinement process should take place at much smaller 
value of T. For T > 100 MeV numerical deviations from the behaviour (|43|) 
are therefore expected, as we will discuss below. 



We can compare these results with those obtained by the authors of [13 
in the framework of chiral perturbation theory at finite temperature within 
the hard thermal loop approximation: 



W(T) = W(0) 



TiaT 



(45) 



A part from the difference in 5m? {Q) discussed above, we get a small devia- 
tion also in the correction oc T^: whereas the term +TTaT'^ is identical, the 
remaining part oc differs by the factor 



3 + /i^ ml + fi^ 
— I In — 

4 4yU2 77j2 



(46) 



whose numerical value is 1.3 at /i = 2 GeV and 1.5 at /i = 3 GeV. This 
difference arises from the different treatment of the virtual photon effect: in 
particular the logarithmic n dependence in ( |43| ) arises from the direct photon 
coupling to pions, which is allowed by the pion electromagnetic form factor 
at T 7^ 0, as given by (PO]). In particular if one assumes that also this direct 

coupling, i.e. the last term proportional to in (|20D, is multiplied by the 

.2 \ -1 



factor ( 1 — ^ 1 , then the same result as in is obtained. 



m 
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Let us now turn to the numerical results of our analysis. In Fig. 1 we 
report two curves. The solid line is the result of the full numerical analysis 
contained in Section ^, while the dashed line represents 5m? {T) as computed 
by the approximate formula (^31). Both curves are obtained at /i = 3 GeV. 
The comparison shows that: 

1) For small T, the small deviation (< 7%) is mainly due to a chiral symmetry 
breaking term which is taken into account in the full calculation (solid line) 
and not in the approximate calculation (dashed line). This extra term is part 
of the pion pole contribution to the dispersion relation (the uj contribution 
is always very tiny). 

2) For large T, i.e. T > 100 MeV, the difference is mainly due to the fact 
that the approximation (^) gets worse with the increasing temperature; in 
particular the positive contribution, i.e. the last term in (|43|), is quadratic 
only for T < 100 MeV, while for T > 100 MeV, it increases more slowly and 
the effect of the negative term remains unbalanced. For T = 200 MeV this 
correction is of the order of 80%. 

All these results are almost independent of in the range 2 — 3 GeV. 

5 Conclusions 

We have extended the Cottingham formula at finite temperature. This for- 
mula allows the computation of time ordered products of two currents be- 
tween hadronic states; even though it is generally applied to the calculation of 
electromagnetic mass differences of hadrons, its possible range of applications 
is wider and includes the evaluation of matrix elements of other products of 
currents such as those occurring in the analysis of weak non leptonic decays. 
Therefore its extension to T 7^ may permit the study of finite tempera- 
ture effects also for weak hadronic processes. As an application, we have 
considered the finite temperature 7r+ — vr*^ electromagnetic mass difference; 
the use of the Cottingham formula leads to results similar to those reached 
by chiral perturbation theory for small temperatures < 100 MeV. We have 
also computed deviations due to chiral symmetry breaking, which are of the 
order of 7% or less, and corrections to the hard thermal loop approximation 
whose role is relevant for T > 100 MeV. 
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Figure caption 

Fig. 1. The difference 5m? = — m^o as a function of tlie temperature. 
Tlie solid line is the result of the full calculation (with a value of the ultra- 
violet cut-off /i = 3 GeV), the dashed line gives the result in the chiral and 
small temperature limit. 
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